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STRUCTURE FUNCTIONS OF MANY-BOSON 
SYSTEM W^ITH REGARD FOR DIRECT THREE- 
PACS 05.30.Jp, 61.20.-p AND FOUR-PARTICLE CORRELATIONS 


On the basis of the expression for the density matrix of interacting Bose particles in the 
coordinate representation with regard for the direct three- and four-particle correlations 
[I.O. Vakarchuk and O.I. Hryhorchak, J. Phys. Stud. 3, 3005 (2009)], the two-, three-, and 
four-particle structure factors of liquid ‘^He in a wide temperature interval were calculated in 
the approximation “one sum over the wave vector”. In the low-temperature limit, the expres¬ 
sion obtained for the two-particle structure factor transforms into the well-known one. In the 
high-temperature limit, the expressions for the two-, three-, and four-particle structure factors 
are reduced to those for the ideal Bose gas. The results obtained can be applied to calculations 
of the thermodynamic functions of liquid ‘^He and to the determination of the temperature 
dependence of the first-sound velocity in a many-boson system. 
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1. Introduction 


The researches of structure functions play an impor¬ 
tant role in studying the Bose and Fermi systems, be¬ 
cause the results obtained theoretically can be com¬ 
pared directly with experimental data. The central 
role in the structural researches of those systems be¬ 
longs to the total scattering cross-section, which is 
called the dynamic structure factor. This parameter 
makes it possible to determine both the spatial struc¬ 
ture of the substance and the structure of its energy 
spectrum [l|,@|. With its help, as well as with the help 
of its derivatives, a lot of different systems are studied 
today, e.g., the Bose gas in a trap [3, linuid "^He (3| 
and ^He y in two dimensions, solid ^He y , thin films 
7|, Lennard-Jones rarefied gas Q, superfluid helium 
sj, parahydrogen 0 nodels with turbulence 0, 
and so forth. 

Besides the dynamic structure factor, not less im¬ 
portant is its zeroth moment or the static struc¬ 
ture factor, which has been measured a lot of times 
in a wide temperature interval. The researches were 
carried out at the saturated vapor pressure within 
the neutron [13 and X-ray (l3| diffraction meth¬ 
ods. Experimental works on the structure factor mea¬ 
surement were analyzed, e.g., in work 0 , where cor- 
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rections were also proposed in order to coordinate 
results of various authors. The Monte-Carlo method 
was applied to study the region around the structure 
factor peak in a Bose condensate, and the peak was 
shown to locate higher than the theoretical value ob¬ 
tained in the framework of the low-density approxi¬ 
mation [l5| . The contribution of three-particle corre¬ 
lations to the structure factor of liquid ^He was found 
in work 0, and a procedure of calculation of the 
effective pair potential on the basis of experimental 
data obtained for the structure factor and with the 
help of the Monte-Carlo simulation scheme was pro¬ 
posed in work (l^ . Among structure functions, we 
also mention the pair correlation function, which is 
one of the key quantities characterizing the coherent 
properties of a Bose condensate 0|- 

In this work, we aimed at finding not only the pair 
structure factor, but also expressions for the three- 
and four-particle structure factors in the approxi¬ 
mation “one sum over the wave vector”. This result 
should help us, in turn, to simplify calculations of 
the thermodynamic functions of a Bose system in the 
approximation “two sums over the wave vector” and 
to facilitate the solution of the still unsolved task to 
describe such a system as liquid “^He in a wide tem¬ 
perature interval and, especially, in a vicinity of the 
A-transition. 
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The structural properties of liquid ^He at low tem¬ 
peratures have been discussed for a long time in the 
framework of the collective variable approach [H- 


23|. However, a theoretical calculation of the pair 


structure factor in a wide temperature interval has 
been carried out only recently in work (^ . by taking 
advantage of the averaging with the density matrix 
of interacting Bose particles. Later on, the structure 
functions in a wide temperature interval were also 
described in other works [ 2 ^. However, the authors 
of the cited works calculated the average values in 
the framework of the density matrix approach and 
in the pair correlation approximation, so that the 
pair structure factor was obtained in the same ap¬ 
proximation. The agreement with exp erimental data 


for the pair structure factor [1^[13,126| is good in 


this case, but incomplete, because, as is known, the 
contribution of many-particle correlations to the ob¬ 
served quantities of a many-boson system can turn 
out rather substantial dilillii. 

In works [ 2 ^, the irreducible two-, three-, and four- 
particle structure factors, as well as the pair distribu¬ 
tion function, were calculated in a wide temperature 
interval making allowance for only the indirect three- 
and four-particle correlations. The obtained theoret¬ 
ical results can be improved by taking direct correla¬ 
tions into account as well. However, in this case, the 
indicated quantities have to be calculated with the 
density matrix containing not only the pair, but also 
three- and four-particle direct correlations. This task 
is a purpose of this work. In our calculations, we will 
base on the approaches proposed in our earlier works 
nils and the results obtained there; in particular. 


these are the expressions for the density matrix and 
the partition function for a many-boson system in a 
wide temperature interval and the methods of their 
calculation in the approximation “two sums over the 
wave vector”. 

An important feature of this work is a graphic pre¬ 
sentation of the results obtained. As a rule, numeri¬ 
cal calculations are carried out for this purpose. The 
input data at such calculations include experimental 
results for the structure factor extrapolated to the ab¬ 
solute zero temperature. The general scheme of spec¬ 
ulations on this topic and the corresponding results 
can be found in work . Continuing the issue of nu¬ 
merical calculations, it is worth paying attention to 
work (^ . where the interatomic interaction poten¬ 
tials were restored on the basis of experimental data 


(as was done in work (^1. and the thermodynamic 
and structural properties of ^He were studied. 

The numerical calculation of the pair structure fac¬ 
tor was carried out taking the effective mass into ac¬ 
count. The expression for the latter was given in work 
( 3 ^ . A necessity of its introduction was substantiated 
in work (^ . 

The expression obtained for the two-particle struc¬ 
ture factor transforms into an already found one 
for the low-temperature limit [^. In the high- 
temperature limit, the expression for the two-, three-, 
and four-particle structure factors are reduced to 
the corresponding structure factors of the ideal Bose 
gas. The two-particle structure factor obtained in the 
approximation “one sum over the wave vector” also 
opens a way to finding the temperature dependence 
of the first-sound velocity in liquid "^He and to com¬ 
paring it with experimental data. 


2. n-Particle Structure 

Factors for a Many-Boson System 

According to the definition, the n-particle structure 
factor equals 

where N is the number of particles, pq = 
— are collective variables, and the 

notation (...) means the averaging with the density 
matrix of interacting Bose particles. In the calcula¬ 
tions to follow, the density matrix in the approxima¬ 
tion “two sums over the wave vector” with the fac¬ 
torized density matrix of the ideal Bose gas will be 
used. This approximation involves the direct three- 
and four-particle correlations, and the matrix itself 
looks like 


R{p\p') = R%{t\t')PMp')P{p\p'). 

where R^(r\r') is the density matrix for noninter¬ 
acting Bose particles, Ppr{p\p') a factor taking into 
account pair correlations, and P{p\p') a factor taking 
the direct three- and four-particle correlations into 
account. In particular. 


Kirlr) = 

X ^exp 

Q 


m 


, 

A^! \2TTph?J 


3N/2 


N 
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where the summation over Q means the summation 
over all permutations of particle coordinates. The fac¬ 
tor which makes allowance for pair correlations looks 
like 0 


Pprip\p') = exp 


where 


bo = -/3Eo + - In 

q^O 


q#0 


'^bi{q)Pq'p~ci- 


q^O 


aq tanh 


tanh 


{<¥) 


1 — e r 


bM = U-: “ 


1 


2 \sinh(/3£'q) sinh(/3£g' 
1 
2 


b 2 {q) = 7, (a, coth(/3£;,) - coth(/3e,)), 


1 / 1 + 


27V 


2m ’ 


X exp 

2 

^7 n 


^0 — ^ '^2, \Pc^P-c^ + Co+2 '^2 C2(q)PqP-q ' 


q^O 




E E E ‘^3*^'^l’^2,q3)PqiPq2Pq3 


qi/0 q 2/0 qa/o 

qi+q2+q3=0 


2 

+ ]v E E C'4(qi,q2)PqiP -qiPq2P-q2 I; (1) 

qi5^0q25^0 

where 

Xq = 262 (g) - bi{q) = 

= aotanh 1 | — tanh (^eq,. 


Explicit expressions for the quantities Cq, (72 (qi), 
C' 3 (qi,q 2 ,q 3 ), and (74(qi,q2) are given in Ap¬ 
pendix 2. They can be obtained using the data quoted 
in Appendix 1. 

3. Pair Structure Factor 

In the case of pair structure factor (n = 2), expression 
o can be rewritten in the form of a derivative with 
respect to the parameter A^: 

(PqP-q) = ^Injy dri... J drNR%ir\r)x 


Eq — EqOq — Ctqj 

vq = j e-*'i'’d>(r)dr 

is the Fourier coefficient of the pair interaction energy 
between particles, and /3 = 1/T is the inverse tem¬ 
perature. An expression for P{p\p') was presented in 
work (^ . Its simplified version can be found in Ap¬ 
pendix 1. As a result, we obtain 

(Pqi-Pq„) = ^/ dri... J dVN X 

X i?^(r|r)Pp^(p|p)P(p|p)pq,...pq„ = 

= ^ J dri... j drAri?^(r|r)pq,...pq„ x 


X exp 


bo - 


q#0 


Piplp)}- 


In the adopted approximation “two sums over the 
wave vector” this expression can be written as fol¬ 
lows: 

(PqP-q) = dri... J drNR%ir\r) x 


X exp 


^0 - q X! ^^q/’-q 


q#0 


— ln{(P(p|p))}. 
dXn 


The expression for the first term was given in works 
(^ . The average {P[p\p)) looks like 

{P{P\P)) = 

_ Jdvi... JdrNR%{ri ,..., rjv|ri,..., rN)Ppr{p\p)P{p\p) 

/dri.../ dr^R^iri, ...,rN\ri, ...,rN)Ppr{p\p) 

It can be obtained on the basis of work (s^ as follows: 

{Pip\p)) = exp|(7o + 2 y]] I ^ + 

I l + A,,bo(gi) 
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>S'o('i'i) 


So{q2) 


qi/O q2#0 


l + Ag^S'o((Zi) l + Aq,2S'o(q2) 


_l_ ^ ^ ^ C*3(qi, q2, q3)-S'o^Hqi, q2, qs) 

n [1 + A5;5'o(gi)] 


qi/O q2/0 qa/o 
qi+q2+q3=0 


i=l 


+fE E 

^'^qi5^0q2 5^0q3#0 «* OWijJJ 

qi+q2+q3=0 

Therefore, using the explicit form for {P{p\p)) and 
the results of works ( 2 ^, we obtain 

S{qi) = (PqiP^qi) = 

^o(gi) 1 


l + A,,^o(gi) [1 + A,,5 o(<Zi)]2 


1 _ C-(4) 


2N 


E 

k2^0 


-^fc2>5'o (qi,-qi,k2,-k2) 


E E 


1 + Afc2iS'o(A:2) 

A/c2 




2-/V 1 + Afe2*S'o(A2) 1 + A/c3S'o((?3) 

qi+k2+k3=0 




. 11 k2, ks 


-4C2(qi)5o'(gi) + 


+ 


C'3(qi) k2, k3)5'Q^^(qi, k2, ka) 
k^OkT^O ^fe2'S'o(fc2)][l + Afc3S'o(fc3)] 

qi+k2+k3=0 

C'|(qi,k2,k3)5o(fc2)5o(fc3) 


^ [l + '''fc2'S'o(fc2)][l + Afe35o(fc3)]y/ 

qi+k2+k3=0 

( 2 ) 


Supposing the terms with a single sum to be small 
in comparison with the quantity corresponding to 
the pair correlation approximation, the two-particle 
structure factor can be written in the form 


Siqi) 


_ >^o(gi) _ jT 

1 + (Agi -f Ilqj )S'o((i'l) 


jjnp 

V 


gi’ 


where 

rrnp ^ ^ (qi.-qi; k 2 ,-k 2 ) 

2NSl{q^)^^ 1 + Afc2^0(fc2) 


2NSl{qr) 


EE 


^/C2 ^k^ 


5('^(qi,k2,k3 


0'^"' k25^0k35^0 
qi-l-k2+k3=0 


[1 + A/c 2 'S'o (^2 )] [1 + Afca (^3 )] 


is the contribution of indirect three- and four-particle 
correlations, and 


n;, = -4G(q.) --Y, 


C4(qi,k2)5o(fc2) 


12 


kSkS 

qi-l'k2+k3=0 


^ k^O 1 + "^fe2'5'o(fc2) 

C'3(qi,k2,k3)5('^(qi,k2,k3) 


[1 + Afe 2 >5'o (^2 )] [1 + Afc3 S'o (^3 )] 


72 C'Kqi, k2, k3)6*o(fc2)<5*o(fc3) 

^ k^Ok^O [l + ^fe25o(fc2)][l + Afc35o(fc3)] 

qi+k2+k3=0 


is the contribution of direct three- and four-particle 
correlations. 

Expression m for the three-particle structure fac¬ 
tor can be presented in the form 


(PqiPqaPqa) ~ 


^/N 6\n{P{p\p)) 

2 (5C'3(qi,q2,q3)' 


A direct calculation on the basis of the previous for¬ 
mula gives the following result: 

-S'^^^(qi,q2,q3) = '/n {pc^^Pq2Pc^3) = 

^ f _ >go^^(qi,q2,q3) _ 

\[1 -I- A53S'o(gi)][l -I- Aq2-S'o((Z2)][l + ^qsSoiqs)} 

_ 12C3(qi, q2, q3)<S'o(gi)5'o(g2)»5'o(g3) | 

[1 -I- Agi5'o((i'i)][l -I- Ag2*S'o((i'2)][l + A53S'o(q3)]/ 


The irreducible four-particle structure factor takes 
the form 


<5'^^^(qi,-qi,q2,-q2) = 

= Af [(pq3P_q,^Pq2P_q2) — (pq^P-qj) (pq2P-q2)]- 


The average (PqiP-qi) was found earlier, and we have 
to calculate (pqjP_qjPq 2 P-q 2 )- Again, on the basis of 
formula ([T]), it can be shown that 


(Pqi P-qi Pq2P-q2 ) 


1 

I\ dXq^ dXq2 


where 



drNR%{r\r) x 
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X exp 


Then, 




q^O 


Pip\p)- 


{PqiP-qiPq^P-q'i) (PqiP-qi) {Pq'iP-qi) ~ 

cP 


dXq^ dXq2 


\iilJdri...JdrNR%ir\r) x 


X exp 


1 


bo - A,PqP_q 


qv^o 


In {Pjplp)) 

dXq^ dXq2 


<S'^'‘^(qi,-qi,q2,-q2) = 


1 


[1 + AgiS'o((7i)]^ 


[1 + Xq^So{q2)]‘^ 


5'o‘‘^(qi>-qi.q2,-q2)- 


2A 


|qi+q2| 


5'o^^(qi,q2,-qi - q2) 


n 2 


1 + ^|qi+q2|'5'o(|qi + q2|) 
48S'o(gi)S'o(g2)5'o^^(qi,q2,-qi - q2) x 
C'3(qi,q2,-qi - q2) 


1 T A|qj+q2|5'o( 

16SUq,)S^iq2) 


qi +q2|) 

C'4(qi,q2) 


18 


C'l(qi,q2,-qi - q2)S'o(|qi + q2|) 


1 + Aiqi+qai'S'odqi + q2|) 


s,(,) = i, = 


dl3 


5'f^(qi.q2,q3) = i, 


^>5'(^^^(qi,q2,q3) 

dp 


= 0 , 
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So (qi,-qi,q2,-q2) = o, 


dp 


= 0 . 


The first term in the expression above was also found 
earlier [2^. The second one is easy to calculate taking 
the explicit expression for {P{p\p)) into account. As 
a result, we obtain 


A straightforward verification also demonstrates that 

lim C 2 (qi) = ia 2 (qi), 

/S-^oo Z 

lim C'3(qi,q2,q3) = ^osCqi, q 2 , qs), 

/3->oo D 

lim C4(qi,q2) = ■^a4(qi,-qi, q 2 ,-q 2 ), 

0^00 O 

where the quantities 02 (qi), a 3 (qi, q 2 7 qs), and 
a 4 (qi) ~qi) q 2 , —q 2 ) are the known expressions [l^ 
and look like 


“ 2 (qi) = ^ E 

q2 5^o 

(q2,qi + q2 


2—a 4 (qi, -qi, q2, -q2) + 


+ 




2gfaqi 

a3(qi,q2,-qi - q2) 


a3(qi,q2,q3) = - 


E (q*qj)(agi - l)(agj “ 1) 

l<i<4<3 


a4(qi,-qi,q2,-q2) = 


2Eq^ 

i=i 

1 


[j^qj 


qfaq^ + qjotq^ 


4. Two-, Three-, and Four-Particle Structure 
Factors in the Low-Temperature Limit 

In the low-temperature limit, the pair and three- 
particle structure factors equal unity, and the irre¬ 
ducible four-particle one equals zero. Their deriva¬ 
tives with respect to the inverse temperature vanish 
in this limit. One may verify it directly by analyzing 
the corresponding expressions. 


X {(qi + q2)^ai(qi + q2, -qi, -q2) + 

+ (qi - q 2 )^a 3 (qi - q2,-qi,q2)- 

- [(qi,q2 + qi)(agi - i) + (q2,qi +q2)(cr92 - i)] x 

X a3(qi +q2,-qi,-q2)- 

- [(qi,qi - q2)(agi - i) + (q2,q2 - qi)(ag2 - i)] x 
xa3(qi -q2,-qi,q2)}- 

Taking the aforesaid into account, we obtain the 
following expression for the pair structure factor in 
the low-temperature limit: 




a 


Qi 


a 


Qi 




^^k2#0k3#0 

qi+k2+k3—0 

k2#0k3^0 

qi-tk2-l-k3=0 


1 Q4(qii k2) 2 a|(qi,k2,k3) 

Iv 2^ TPi. ^Iv 2^ 2-^ 


N ^ Ofc 
k2/0 


N CXkoCtk-i 

k2#0k3^0 

qi+k2-l-k3=0 


, (3) 


1119 


































I.O. Vakarchuk, O.I. Hryhorchak 


where a 4 (qi,k 2 ) is an abbreviated notation for the 
quantity a 4 (qi, —qi, k 2 , —k 2 ). In the adopted ap¬ 
proximation “two sums over the wave vector” the 
structure factor can be written in the form similar 
to that in work (^ . 


lim C'4(qi,q2) = hm C 4 (qi,q 2 ) = 
p^O p^O o 

so that 

Im (^2(qi) = ^1^(73(qi,q2,q3) = ^1^(74(qi,q2) = 0 


S{q) = l/[l-2a2{q)-n<l)l 

where 

^ 1 ^a4(k,-k,q,-q) 2 


xEi: 

ki5^0 k2#0 
q-|-ki-|-k 2=0 


{d2(fci)d2(fc2)+a3(q, ki, k2)[l-fa3(q, ki, k2)]} 


(1 - 2 d 2 (fci))(l - 202(^2)) 


02 ( 91 ) = - 1) + 


N 


E 

q2/0 


92 


29fa 


-a4(qi,-qi,q2,-q2) 




(q2,qi +q2) 


9?a5i 


a 3 (qi,q 2 ,-qi - q 2 ) 


Analogously, we obtain the following expressions 
for the three- and four-particle structure factors in 
the low-temperature limit: 


1 


5'^^Hqi>q2,q3) = 

^93 

<S'^'‘^(qi,-qi,q2,-q2) = 

4 


{1 -f 2a3(qi,q2,q3)}, 


“^91 “^92 L 


‘^|qi-i-q2| 


- 1 


'^|qi-l-q2| 


-f 


-03(qi,q2,-qi - q2) 


-a3(qi,q2,-qi - q2! 


®|qi-l-q2| 

-f a4(qi,q2) + 

‘^|qi-i-q2| 

5. Two-, Three-, and Four-Particle Structure 
Factors in the High-Temperature Limit 


Therefore, in the high-temperature limit, the two-, 
three-, and four-particle structure factors for a many- 
boson system transform into the corresponding ex¬ 
pressions for the ideal Bose gas: 

lim S'( 9 ) = S'o(9), 

im 5'^^Hqi>q2,q3) = -S'o^^(qi,q2,q3), 

1 ™'S'('‘^(qi,-qi,q2,-q2) = 5'^'‘^(qi,-qi,q2,-q2)- 


6. Numerical Calculations 

The numerical calculation of the two-particle struc¬ 
ture factor ([2]) will be carried out taking the effec¬ 
tive mass into account [s^. In order to not exceed 
the calculation accuracy, the effective mass will be 
used only in the terms that reproduce the pair cor¬ 
relation approximation. At the same time, the ex¬ 
pressions containing the sum over the wave vector 
will contain a “bare” mass. Here, the following re¬ 
mark is worth making: in the structure factors of 
the ideal Bose gas which enter the expressions with 
a sum over the wave vector, the effective mass is 
used only to shift the critical point owing to the ac¬ 
tivity renormalization zq = exp[/3p], where p is the 
chemical potential. The introduction of effective mass 
makes it possible to avoid infra-red divergences in the 
non-renormalized four-particle structure factor of the 
ideal Bose gas. 

To calculate the quantities with a single sum over 
the wave vector, we should change from summation 
to integration according to the well-known rule 


Using the explicit expressions for the quanti¬ 
ties C' 2 (qi), C 2 (qi), C'3(qi,q2,q3), C' 3 (qi, q 2 , qa), 
C 4 (qi,q 2 ), and C' 4 (qi,q 2 ) (see Appendix 2), we 
can easily obtain that, in the high-temperature limit 
(T —>• 00 or /3 —>• 0), 


hrn^ C 2 (qi) = ^im ^2 (qi) 


1 


EE 


9?(q2q3) 


I 6 N ^ ^ OoOo 

q2/0q3/0 ^2^3 

qi-l-q2-fq3=0 


lim C'3(qi,q2,q3) = lim ^“(qi, qa, qa) = 

p—p—J. ^ 


E 

k 


V 


(27t)3 


dk. 


After the corresponding transformations and the re¬ 
quired changes in variables, we obtain the following 
rule for the change from summation to integration in 
our case: 


1 

N 


EE 

k2#0 k35^0 
qi-l-k2+k3=0 


1 


191 +^ 2 ! 


47r^pqi 


k^dko 


ksdks, 


\qi-k2\ 
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where p is the equilibrium density in the Bose sys¬ 
tem. For such quantum liquid as ^He, the latter pa¬ 
rameter equals p = 0.02185 A [36|. The next step 
consists in calculating the quantities Uq with the use 
of a pair structure factor at the zero temperature ex¬ 
trapolated on the basis of experimental data. The cor¬ 
responding information is taken from work . 

Now, let us rewrite Eq. ([3]) in the form 




OLq 


1 

^qi L 


-EE 




^k2 


2N 

k2/0 ks/O 
qi+k2+k:3—0 

■2a2(qi) + E E 

k2/0 ka/O 

qi+k2+k3—0 

1 a4(qi,k2) 2 




«3(qi,k2,k3) 
C^k2 ^^3 


k2#0 


^k2 


N 


EE 


o§(qi,k2,k3)' 


k2/0 k3 5^0 
qi+k2-(-k3=0 


Q^/C2 


(4) 


This is an iterative equation for Uq. In the zero-order 
approximation, we have Oq = l/S'®^P(g). Substituting 
this Qfg-value into the right-hand side of equality ([4]), 
we obtain the a^-value in the first approximation, and 
so forth. However, this iteration process does not con¬ 
verge, which is most likely connected with an insuf¬ 
ficient number of terms in the series expansion for 
the structure factor ([3)). Therefore, the consideration 
will be confined only to the zero-order approximation 
for Uq. 

The results of numerical calculations for tempera¬ 
tures of 1.0, 1.38, 1.67, 2.2, 2.5, 3.0, 3.5, and 4.24 K 
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Fig. 5. The same as in Fig. 1, but at T = 2.5 K 




Fig. 6. The same as in Fig. 1, but at T = 3.0 K 



are exhibited in Figs. 1 to 8, respectively. Experimen¬ 
tal data for the structure factor at those temperatures 
were taken from works [l^[T^ . In the presented fig¬ 
ures, the solid curves correspond to the structure fac¬ 
tor calculated with regard for the direct three- and 
four-particle correlations, the dashed curves to the 
pair correlation approximation, and the circles to the 
experimental structure factor values. 

7. Conclusions 

In this work, expressions for the two-, three-, and 
four-particle structure factors in a wide temperature 
interval were found in the approximation “one sum 
over the wave vector” with regard for the direct three- 
and four-particle correlations. In the low-temperature 
limit, the expression obtained for the two-particle 
structure factor transforms into the well-known one 
[^ . The same is valid for the high-temperature limit. 

The derived expressions are rather cumbersome. 
They were analyzed, by using numerical methods, 
and graphic representations of the pair structure fac¬ 
tor at various temperatures of liquid ^He were plot¬ 
ted. The calculation of the internal energy and the 
determination of the temperature dependence of the 
first-sound velocity in the many-boson system will be 
a subject of our next papers. 


APPENDIX 1 

P(pIp') = exp 


+ E E 

qiT^O ii=0 ji=0 
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+ ^EE E E c4(P3,-r3,2^^-2^3); 
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2sh 
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The superscripts 2i, ^2, 23, Ji, and j2 at the quantities 
pq can acquire values of 0 or 1, namely, 0 means the 
absence of the prime, and 1 its presence. The notation 

C2 —1^1) stands for C2 (^pgi, , 03(1*1,2*2,3^3) 

for C3(pg\,pg,pg), and C4(l-?'i,-l*i, 2-?'2 ,-2*2) for 
Ci{pn , , P92> P-q^)- Accordingly, 

C2 (l-’\-lA) =C2 -c§ (1^'1,-lA), 

C3 (iA, 2*2,3'3 ) = c 3 (iA, 2*2,3*3) -cg(lA,2*2,3*3), 

C4(Pi, -lA , 2^2 ^ _ 2 i 2 ) = c 4 (iA ,-r\ 2 ^^ , -2*2) - 
-c^(iA,-iA^2'’2,-2*2). 

Here, the quantities Cq, c®, C3, and mean cq, C2, C3, and 
C4, respectively, in which aq^ = cnq^ = 0:^3 = 1. The quantity 
C4(l-^i, —1*1,2-^2 ^ —2*2) consists of two terms, 

C4(l-^‘i, -1^1,2-^'2 , -2*2) = (1^1, -1*1,2-^2 ^ _2^2 ) + 
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H--—— - Q{—aq^ , —092 > “ 93 )^ 

EyEq-i + tjq^) 

X ch — (2(12 — l)Eq2 + 2 (il — l)Eq^ — Eq^ J-. 

The coefficients C2(l'’i, —I'l) and cq can be ex¬ 
pressed in terms of (fn , —l®i, 212, —2*2) and 

(111,-1 A, 212,-2*2) as follows: 

C2(l'’^ , —I'l) = 20:92 [ 0 rq 2 ] X 

X ( 2 c[^'> (fii, -1*1,2, -2') -I- cf^ (111, -1*1 ^ 2, -2')), 
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The quantity 03(1*1,2*2,3*^) looks like 

C3(1*1 2*2 3*3)- ^ (-i)^^+^"(qiq2) ,, 

’ ’ 16 2m sh[/3Eg^] shf/SEgg] shf/SEgg] 


^EE (Q 91 Q 92 ±1 ± 2 !) 
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[2 -^92 + ^Q3 
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X sh ((-1)*1E91 + (-1)*"E92 + (-1)’==-E93jj • 

In the expressions written above, the following notations were 
introduced; 

Eqi = ltli?93; Eq2 = lt2-f^91 i 

d93=±109j; Q92=±ia92; 

E = ±1-E93 ±2 Eq 2 + Eq^', 

(3(091,092,093) = (±1 ±2 O91Q92 -b l)(qiq2)-|- 
-I- (±1091093 ± l)(qiq3) ± (±2092093 ± I)(q2q3). 
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APPENDIX 2 

C 2 (qi) = C 2 (qi) - C 2 (qi), 

C3(qi,q2,q3) = C3(qi,q2,q3) - C 3 (qi, q2, qs), 
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1.0. BaKapnyK, O.I. FpmopnaK 

CTPYKTyPHI OyHKIl,!! 
BAFATOBOSOHHOI CHCTEMH 
13 BPAXyBAHHHM HPHMHX TPH- 
TA HOTHPHHACTHHKOBHX KOPEJIHUIH 

P e 3 K) M e 


Ha ocHoni Bupasy MaxpHu;! rycxHun BsaeMO^iiouHX 6o3e- 
uacxHHOK B KOop^HuaxHOMy 3o6paxceHHi i3 BpaxyBaunaM npn- 
MHX xpu- i HOxupHuacxHHKOBHx Kopejmu;iH [1. O. BaxapuyK, 
O. 1. PpHropuaK, HCypu. 4>i3. ;],ocji. 3, 3005 (2009)] 6yjm po3- 
paxoBaui abo-, xpn- i uoxHpHuacxHHKOBi cxpyKxypui 4>aKxo- 
pH pi^Koro “^He b HabjinxceHni “o^mei" cyain 3a xbhjibobhm 
BeKxopoM” ^jm mupoKoro iuxepBajiy xeMnepaxyp. B rpaHHu;i 
HH3BKHX xeMnepaxyp oxpHMauHH Bupaa AJin ;],BOHacxHHKOBoro 
cxpyKxypuoro 4>aKxopa nepexo;],Hxi> b yxce Bi^OMnii. B rpaHHu;i 
BHCOKHX xeMnepaxyp BnpasH A-nn abo-, xpn- i HOXHpnuacxHH- 
KOBHX cxpyKxypnnx 4)aKxopiB peAynyioxBCn ao cxpyKxypnnx 
(^aKxopiB iAoajiBHoro 6o3e-ra3y. PeayjiBxaxn poSoxn Mo:>KyxB 
6yxH aacxocoBani a^i^ poapaxynxy xepMOAHHaMiuHHx 4^yHKTT,in 
piAKoro “^He i 3HaxoA2*ceHHn xeMnepaxypnoi’ aajiexcHOCxi mBHA- 
Kocxi nepmoro 3ByKy b 6araxo6o30HHiH cncxeMi. 
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